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. We perform a study of the flavour evolution in the early universe of a multi-flavour active-sterile 

^ ' neutrino system with parameters inspired by the short-baseline neutrino anomalies. In a neutrino- 

symmetric bath a "thermal" population of the sterile state would quickly grow, but in the presence of 
, primordial neutrino asymmetries a self-suppression as well as a resonant sterile neutrino production 

can take place, depending on temperature and chosen parameters. In order to characterize these 
effects, we go beyond the usual average momentum and single mixing approximations and consider 
a multi-momentum and multi-flavour treatment of the kinetic equations. We find that the enhance- 
ment obtained in this case with respect to the average momentum approximation is significant, up 
' to ~ 20% of a degree of freedom. Such detailed and computationally demanding treatment further 

raises the asymmetry values required to significantly suppress the sterile neutrino production, up 
to \L^\ > O(10~^). For such asymmetries, however, the active-sterile fiavour conversions happen 
(~| ' so late that significant distortions are produced in the electron (anti)neutrino spectra. The larger 
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increasing the *He abundance in primordial nucleosynthesis (BBN). The standard expression of 
the primordial yields in terms of the effective number of neutrinos and asymmetries is also greatly 
altered. We numerically estimate the magnitude of such effects for a few representative cases and 
comment on the implications for current cosmological measurements. 



PACS numbers: 14.60. St, 14.60.Pq, 98.80.-k 26.35.-l-c 
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O ■ I. INTRODUCTION 

o 

In recent years different short-baseline neutrino oscillation experiments have found anomalous results that may 
^ be interpreted by enlarging the standard description of neutrino oscillations in terms of three active species. In 
^ particular, the — oscillations in LSND [Ij and MiniBoone 2] experiments (recently constrained by the ICARUS 
. experiment Q), the Ve and i^e disappearance revealed by the Reactor Anomaly 141, and the Gallium Anomaly Q 
I ' can be described in terms of light [m ~ C(l)eV] sterile neutrinos which mix with the active ones (se e 0, Q for 
ILJ , recent reviews). In this context, scenarios with one (dubbed "3-1-1") or two ("3-f2") sterile neutrinos [8l-[l4| have been 
. ^ ■ proposed to fit the different data. 

] Cosmological measurements represent a powerful tool to probe the number of neutrinos and their mass at cV 
?H ■ scale (see, e.g., [lB4I3)- The non-electromagnetic radiation content in the universe is usually expressed in terms of 
the effective number of excited neutrino species, NeS. This can be constrained by Cosmic Microwave Background 
(CMB) il8H20], Large Scale Structure (LSS) ^^^^ and Big Bang Nucleosynthesis (BBN) data [UHll. Current 
measurements, especially CMB ones, slightly favor the existence of some extra radiation, though the results of WMAP- 
9 data release [l^ and Atacama Cosmology Telescope [2?'], when combined with other cosmological constraints, are 
compatible with the Standard Model expectation value, A'off = 3.046 [28|. In particular, WMAP-9 finds A'off = 
3.84 ± 0.40 when the full data are analyzed [2^, while Atacama gives Ncs — 2.79 ± 0.56. In this context, a recent 
Bayesian analysis of the current cosmological datasets does not show a strong preference for a value of Ncs larger 
than the standard one [25j . Very recently, a combination of cosmological data notably including the ones released by 
Planck has yielded A^off = 3.30 ± 0.27 

Similarly, up to about one fully thermalized sterile neutrino is still marginally allowed by BBN [2l|, , while a 
value corresponding to two sterile states appears largely excluded [1^, but no significant preference for a larger- 
than-standard value of A^cff either. Notice, however, that even a single extra thermalized sterile neutrino with mass 
m ~ 0(1) eV appears to be inconsistent with mass bounds from CMB and LSS data [H, Isol - fssl lssj . 

In order to reconcile the eV sterile neutrino interpretation of the short-baseline anomalies with the cosmological 
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observations, the most straightforward possibility would be to suppress the sterile neutrino thermalization in the early 
universe, correspondingly reducing the expected excess in extra-radiation. In this sense, it was proposed at first in 
in [sBl (see also [s^l) to consider a primordial asymmetry between neutrinos and antineutrinos ^ 

= ^ . (1) 

In principle, one would expect the neutrino asymmetry to be of the same order of magnitude of the baryonic one, 
•q = (riB — njj)ln-y ~ 6 x 10^^". This of course holds for charged leptons, due to the stringent requirement of 
charge neutrality of the universe, but not necessarily for neutrinos. In fact, the constrains on are quite loose, 
allowing also \L^\ ~ 10^^ — 10~^ (solfisj . A primordial neutrino asymmetry would add an additional "matter term 
potential" in the active-sterile neutrino equations of motion. If sufficiently large, one expects this term to block the 
active-sterile fiavour conversions via the in-medium suppression of the mixing angle. Nevertheless, this term can 
also generate Mikheev-Smirnov-Wolfenstein (MSW)-like resonant flavour conversions among active and sterile 
neutrinos, enhancing their production. In order to assess which of the two effects dominates the flavour evolution it 
is mandatory to perform a study of the kinetic equations for active-sterile neutrino oscillations. 

In a recent paper [50l | some of us performed a flrst study of active-sterile flavour conversions in the presence of 
neutrino asymmetries, considering (3-1-1) and (2-1-1) scenarios inspired by the recent flts of all the short-baseline 
neutrino anomalies. In order to simplify the numerical complexity of the problem, we adopted equations of motion 
integrated over momenta, often referred to as "average (or single) momentum approximation" . Loosely speaking, 
this can be thought of as an approximation in which all neutrinos share the common thermal average comoving 
momentum (p)/T = 3.15. Under this assumption, we found that in the case of equal asymmetries among the active 
species, a value of \Li,\ ~ 10"'^ is required to start suppressing the resonant sterile production. An even larger results, 
~ 10~^, is necessary to lower the sterile neutrino abundance in case of opposite initial asymmetries. 

However, A^eff is not the only parameter that can affect the cosmological observables. Indeed, for such large 
values of resonant active-sterile neutrino conversions occur near or after the decoupling temperature for the active 
neutrinos, making ineffective the repopulation of active species through collisions. The lack of repopulation of electron 
neutrinos would in general produce distorted distributions which can move up the n/p freeze-out and hence increase 
the ^He yield, the main product of BBN. In order to characterize the possible distortions in the active neutrino 
spectra, it is necessary to go beyond the average momentum approximation and consider a detailed treatment of the 
full momentum-dependent kinetic equations, due to the momentum-dependence of the resonant conversions between 
active and sterile neutrinos. 

The purpose of the present work is to perform for the first time a full multi-flavour and multi-momentum treatment 
of the active-sterile neutrino oscillations for the (2-1-1) scenario considered in [50l| in the presence of primordial neutrino 
asymmetries. It was shown in 50] that this model captures the main features of the complete (3-f 1) scenario. In 
particular, we follow the evolution of the system in the presence of \Lv\ ^ 10~^ where the distortions of the active 
neutrino spectra start to become sizable. We remark that for smaller values of the active neutrino asymmetries 
the thermalization of the sterile neutrinos is complete, producing a tension with cosmological data. In the studied 
cases, we observed an enhancement in the sterile neutrino production of up to AA^cg' ~ 0.2 with respect to what 
observed in the average momentum study [s^. This implies that one needs to consider even larger asymmetries 
(at least of the order of \L^\ ^ in order to significantly suppress the production of sterile neutrinos. On the 

other hand, for these values of the asymmetries we find relevant distortions in the electron (anti)neutrino spectra 
notably modifying the BBN yields (see, e.g., [Sll, [s^l)- Computing reliably these distortions and A"eff as functions 
of the asymmetry parameters is a very challenging task, involving time consuming numerical calculations for the 
fiavour evolution. Nonetheless, the few representative values of the asymmetries we investigated already allow to 
draw non-trivial implications for cosmological observables. 

This article is structured as follows: in Section II we briefly summarize our formalism, highlighting the modifications 
with respect to [soj . We present our results concerning the sterile neutrino production in the multi-momentum scenario 
and the differences with respect to the average momentum case. Then, in Section IlIII we show the distortion of the 
electron neutrino spectra for different values of asymmetries and we discuss the impact on observables, namely Aoff 
and ^He and abundances. Finally, in Section ITVl we conclude. 



Alternative escape routes have been recently proposed, see e.g. [si 
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II. MULTI-MOMENTUM FLAVOUR EVOLUTION 
A. Equations of motion 

In this Section we summarize the Equations of Motion (EoMs) for the (2+1) active-sterile neutrino system in the 
early universe, using the same notation of 50], to which we address the reader for details. In particular, we describe the 
time evolution of the neutrino ensemble in terms of the following dimensionless variables (replacing time, momentum 
and photon temperature, respectively) 



ma 



y = pa 



(2) 



where m is an arbitrary mass scale which we put equal to the electron mass nif,. Note that the function a is normalized, 
without loss of generality, so that a(t) — >■ 1/T at large temperatures, T being the common temperature of the particles 
in equilibrium far from any entropy-release process. With this choice, can be identified with the initial temperature 
of thermal, active neutrinos. 

In order to take into account the interplay between oscillations and collisions of neutrinos, it is necessary to describe 
the neutrino (antineutrino) system in terms of 3 x 3 density matrices g (g)"^ 



gee ge^i ^?es 

q{x., y) = I g,,e gfj^f, gfj.s 

Qse Qsfi Qss 



(3) 



In this formalism we write the EoMs for g and g as [39|, |53|, |54 
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(5) 



which has to be solved along with the covariant conservation of the total stress-energy tensor 



de „^ 
X— ^e-3V. 
dx 



(6) 



In the previous expressions H denotes the properly normalized Hubble parameter 



H = —H = 



x"^ 8ne{x,z{x)) 
m Y 3M|,; 



M 



PI 



m \ 8Tre{x,z{x)) 



(7) 



and the total energy density and pressure enter through their "comoving transformed" values e = e(x/m)^ ~ e^+e^+ei, 
and V = P{x/m)'^. Compared with the treatment in (50| . we now take into account the non-relativistic transition 
of the electrons and the entropy-transfer to the photons, responsible for increasing the temperature of photons with 
respect to the one of neutrinos, i.e. for making z{x) > 1. However, since we are only interested to situations where 
the change with respect to the instantaneous decoupling limit value I%s = 3 is much larger than the ~ 1.5% effect 
due to incomplete neutrino decoupling during — e~ annihilation |28j . we neglect the latter effect. We can thus 
track z{x) once and for all by using the integrated entropy ratio formula (see e.g. Eq. (15) in |:55], with dfi/dx at the 
right-hand side put to zero). 

The first term at the right-hand side of Eqs. (HI)-® accounts for vacuum oscillations, where in the flavour basis 



= WM^U. UereU =U( 



' e^ ; *^es 1 ^ 



) is the 3x3 active-sterile mixing matrix, parametrized as in [50| . We assume 



Here refers generically to a non-electron active flavour state. 
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9e^i equal to the active 1 — 3 mixing angle Q12, [56|, while we fix the active-sterile mixing angles to the best- fit values 
of the different anomalies [l3] 

sin^^e^ - 0.024 , (8) 
sin2 6'e, = 0.025 , (9) 
sm^Of^s = 0.023 . (10) 

The mass-squared matrix M.^ — diag(— Am^t,r, /2. -I- Am^^^ /2, ArngJ is parametrized in terms of the atmospheric 
mass-squared difference Am^^^ = 2.43 x 10~^ eV^ (S^ and of the active-sterile mass splitting Am^j = 0.89 eV^, fixed 
from the short-baseline fit in 3-1-1 model [l^. In the following we assume the normal mass hierarchy Am^tjjj > 0. We 
checked that results similar to the ones we will present would have been obtained considering the (2-1-1) sub-sector 
associated with the solar mass squared difference Am^^j and with the 1 — 2 mixing angle 9i2. 

The terms proportional to Gi? in Eqs. ©-([S]) encode the matter effects in the neutrino oscillations. In particular, 
the term is related to the energy density of electrons and positrons (see Eq. (20) in ^5^]), while the v — v interaction 
terms are given by 

N,. = ^ j dyy^{Gs{Q{x,y)-Q{x,y))Gs + '^sT:T:[{Q{x,y)~Q{x,y))Gs\}, (11) 
Ez. = ^ y" dyy'^Gs{g{x,y) + Q[x,y))Qs. (12) 

Note that the matrix is related to the difference of the density matrices of neutrinos and antineutrinos, while Ej^ 
is related to their sum. The matrix = diag(l, 1,0) in flavour space contains the dimensionless coupling constants. 
We remark that considering more than one active species, the matrix also contains off-diagonal terms. In the 
presence of large neutrino asymmetries among the active species, as the ones we are considering in this work, the 
matter effects are dominated by the contribution. This latter term makes the EoMs non- linear and are the main 
numerical challenge in dealing with this physical system. At large temperatures it dominates over the vacuum term, 
suppressing any flavour conversions. However, as the universe expands a resonance condition can be satisfied at lower 
temperatures, as extensively illustrated in [50| . 

Finally, the last term in Eqs. (|4])-([5]) is the collisional one proportional to for which we adopt the same 
approximate expression C[p] as in Eq. (28) of [s^l, but now keeping the y— dependence, see [111. This guarantees 
that: i) the correct collisional term is recovered when integrating the EoMs for Q{y) over momenta, ii) the overall lepton 
number conservation is preserved. Note that this is not the case for alternative damping prescriptions often found 
in the literature, where the lepton number conservation is achieved imposing an additional equation (see e.g. 59])- 
We remark that possible minor inaccuracies in the y— dependence of the collisional terms are of little concern for our 
application, since in the cases where large spectral distortions are produced in the electron (anti)neutrino distributions, 
these are caused by oscillatory terms rather than collisional ones. In other words, whenever collisional repopulation of 
depleted active neutrinos is relevant, one has very little departures from thermal spectra, while in the opposite limit, 
the major distortion is not due to the exact y— dependence of the collisional terms. On the other hand, implementing 
an operator enforcing lepton number conservation (numerically to a very high degree of precision) is crucial to avoid 
spurious numerical effects. 

In order to fix the initial conditions for the flavour evolution, we notice that active neutrinos are produced in 
the very early universe with their energy spectrum kept in chemical and kinetic equilibrium by weak interactions 
until temperatures T ~ few MeV, when the corresponding collision rates fall below the cosmological expansion rate. 
In the presence of primordial neutrino asymmetries, the original active neutrino spectra are given by Fermi-Dirac 
distributions parametrized in terms of a temperature Ti, and chemical potentials /i^ for a — e,fj,,T. Each neutrino 
asymmetry in Eq. ([I} can be expressed in terms of the corresponding degeneracy parameter = ^J■a/Tl, as 

with C(3) ~ 1.202, where the right hand side numerical expression corresponds to the leading order in small ^q.. In 
the rest of the paper, we shall indicate the neutrino asymmetries in terms of the parameters rather than L^, in 
order to conform with the more frequently used notation in phenomenological papers. 
The initial conditions for the density matrix g are then given by 

Qin = diag(/oq(?;,Ce),/cq(?;,^M)'0) ' 

Pin = diag(/oq(y,-Ce),/oq(y,-CM)'0) ' (14) 
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with /eq(y, e) = l/[exp(y - e) + 1]. 

We remark that since in our study we consider initial distributions for active neutrinos cfose to their equihbrium 
ones ^, the oscillations among the three active species have a qualitatively important role but quantitatively their 
details are of minor relevance for the evolution of the sterile neutrinos. Therefore, the (2+1) scenario we consider is 
a good proxy for the complete (3+1) scenario. We also mention that in a recent work [571 ) . it has been performed a 
multi-momentum study of the kinetic equations for a system consisting of only one active and one sterile neutrino 
species [i.e. (1+1) scenario]. While instructive in several respects, this is quite a simplified scenario since it does not 
allow to incorporate effects like the existence of more than one active-sterile mixing angle and different choices of the 
neutrino asymmetries for the different flavours. 



B. Results 



In this Section we present our results for the sterile neutrino abundance in the (2+1) scenario described in Sec. Ill Al 
We numerically solve the EoMs Eqs. (2])- ([5]) with an integration routine for stiff ordinary differential equations taken 
from the NAG libraries [6l| and based on an adaptive method. The range for x is chosen to be a; € [2 x 10~'^,0.5]. 
As a compromise between energy resolution of the spectral distortions and computational cost, we took Ny — 21 
momentum modes in the range y S [0, 10]. The grid points are not chosen to be equally spaced, but are instead fixed 
by imposing weighted Gaussian quadrature of the integrals in the right-hand-side of Eqs. (|lip - P^ . By increasing 
the momentum grid points to A^j^ = 30, we checked in some test runs that this is enough to keep the error below the 
per-cent level on the effective number of neutrinos Nee- 

We remark that due to the momentum dependence of the resonance conditions, in the multi-momentum treatment 
of the EoMs there can be significant deviations with respect to the evolution predicted by the average momentum 
scheme. A direct comparison between the single-momentum and the multi-momentum results is reported in Fig. [TJ 
which shows: (a) the momentum- integrated sterile neutrino density matrix element (solid curves) normalized to the 
integral of a Fermi-Dirac distribution with zero chemical potential, 

Jdy y^Qss{x,y) 
'-^^^ ^ Jdy y^UM ' 

(h) the sterile neutrino density matrix element p^g in the average momentum scheme (dot-dashed curves), normalized 
correspondingly. In the left panels we take equal initial neutrino asymmetries, = C/j, while in the right panels we refer 
to opposite ones, = In the upper panels we consider = 10"'^, while in the lower panels we take = 10~^. In 
all the considered cases the values of pss for the single momentum underestimates the sterile neutrino abundance with 
respect to the multi-momentum case. The enhancement obtained with the multi-momentum treatment is significant, 
roughly ~ 20% of a degree of freedom. Moreover, we note that the sterile production in the multi-momentum case 
occurs at higher temperatures with respect to the average momentum case. This is due to the fact that in the 
multi-momentum evolution the sterile neutrino population can start building up earlier via lower momenta modes, 
that resonate earlier than the average momentum. This anticipates the dynamical evolution of ^ and is the main 
responsible for the difference observed. Also, we confirm the observation of [62j that the resonance is more adiabatic 
at higher temperature. Hence, the average momentum treatment of the EoMs is generically expected to underestimate 
the sterile neutrino abundance. 



III. IMPACT ON OBSERVABLES: N^s AND LIGHT NUCLEI ABUNDANCES 



Phenomenological quantities affected by active-sterile neutrino fiavour conversions notably depend on the overall 
non electromagnetic radiation content, parametrized via iVoff , and from the distortions of the electron (anti)neutrino 
spectra, the latter being a basic input for BBN weak rates ^. 

In the upper panels of Figs. 2 and 3 we show the j/-dependent Ve energy spectrum y'^Qeeiy) (dashed curve) at 
T = 1 McV, compared with the initial one y^/oq(y,Ce) (solid line). In particular Fig. 2 refers to = lO^'^ and Fig. 3 



^ This is conservative since some population of sterile neutrinos cannot be excluded in extensions of the Standard Model, following for 

example from decays of heavier neutrino singlet states. 
* Actually, treating properly effects sensitive to the neutrino masses require knowing the flavour composition of the neutrino ensemble as 
well. Similarly, precision computations of CMB anisotropics are in principle sensitive to the neutrino phase-space distributions, see e.g. 
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FIG. 1: Evolution of the total value of the sterile neutrino density matrix element pss as function of the temperature T for 
the multi-momentum case (continuos curves) and the average momentum (dot-dashed curves) case with a thermal momentum 
(y) = 3.15. Left panels correspond to ~ (./i, while in the lowers = —Cm- Upper panels refer to = 10"'^ and right to 



to = 10 ^. In each Figure, in the left panels = Cm' '^hile in the right ones = — C/i- In order to characterize the 
distortion in the spectra with respect to the initial one, in the lower panel we plot the ratio 

Qee{y) 
/oq(2/,4e) 

In the case of = 10^"^, R ^ 0.95 for equal asymmetries and R ^ 0.98 in the case of opposite asymmetries. Conversely 
the spectral distortions are more evident for — 10^^. Namely, for equal asymmetries one finds R ^ 0.82, while for 
opposite asymmetries R ^ 0.9. Indeed, spectral distortions in the active neutrinos are more evident when resonant 
active-sterile conversions occur near the active neutrino decoupling temperature, as pointed out at first in [63j . 

We remind that, as already commented in fs^lj the dynamical evolution of the asymmetries is such that both 
neutrinos and antineutrinos get populated resonantly, roughly in equal values. Hence, we only show here the results 
referring to neutrinos. Of course, in the numerical computation the small differences between the two sectors have 
been properly accounted for. 

Concerning the effective number of neutrino species, we remark that at the level of approximation we are adopting, 
iVcff = 3-1- ANcs enters the dynamics only via its contribution to the Hubble parameter, see Eq. (O, where it rescales 
the standard neutrino energy density contribution as 

3) 1-f — ^ . (17) 



Technically, we compute AN^g in the above equation via the following (numerical) integral 

AN,s ^Y^Jdy y'Tr[^.(.T, y) + g{x, y)] - 2 , (18) 

the factor "—2" being due to the fact that we are considering only two active neutrino species. 

The quantity defined in Eq. (1181) is shown in Figure 2] for two representative values of asymmetries, (^e = 10~^ for 
solid curves and = 10~^ for dashed curves), taken equal (opposite) for the e and /i sector in the left (right) panel. 
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FIG. 2: Cases with = 10 Upper panels: Final Ue energy spectra at T = 1 MeV (dashed curve) and initial ones (continuous 
curve). Lower panel: Ratio R between final and initial energy spectra. Left panels refer to = C/j while right panels are 
for = -Cm- 



For = 10~^, we see that the resonant production of sterile neutrinos starts around T ~ 5 MeV for equal asymmetries 
and T ~ 8 MeV for opposite ones. Since active-sterile neutrino conversions mostly occur when the collisional regime 
is still operative, the active neutrino species are almost fully repopulated, reflecting in a AA'^cff ~ pss in both cases 
(see Fig. 1). This is still true to some extent for the case = — C/x = 10~^, but not when = Cm = 10~^. Namely, 
in this latter situation, the resonant population starts at temperatures as low as T ~ 2 MeV, comparable to neutrino 
decoupling temperature. This means that active neutrino repopulation is only partial. Also, an appreciable difference 
(in this case of ~ 0.1) is established between and N^e. These effects were already noted in [50|, where they were 
even more prominent due to the "less effective" sterile production associated to the single momentum approximation. 

The numerical values of the AiVoff's found in a few representative runs are reported in Table |I] along with the 
values of the yields of "^He mass fraction Yp and deuterium ^H, as obtained from a modified version of the numerical 
code PArthENoPE [6^ for a baryon fraction ujb = 0.02249 and the neutron lifetime r„ = 880.1s, following Particle 
Data Group 2012 recommendations 67,]. For comparison we also show cases with neutrino asymmetries but no sterile 
neutrinos as well as the standard BBN case. 

Technically, note that the rates r„_^p[/,y^, and rp_j.„ , ] are functionals of the distributions If 
we denote with F" the rates computed in the Born approximation for Fermi-Dirac spectra, and with F the actual 
rates for the cases at hand, we have computed the effect of sterile neutrinos by rescaling the rates implemented in 
the code PArthENoPE [66] (see also [g^) by F/F°, which has been numerically evaluated and then interpolated. This 
amounts to a first-order correction in a perturbative approach. However, since the rate corrections is at most C(3%) 
for the largest asymmetries, the error due to this approximation is safely below 0.3% (see e.g. [1^ for the analysis of 
corrections to Born weak rates). Thus, this is comparable or lower than neglecting the modification to the reheating in 
the standard scenario. In determining the shape of the distribution function, especially for the cases with |C| — 10"^ 
the main source of error comes in fact from the discretization of the neutrino distribution (and the corresponding 
interpolation). Still, test runs with a Ny = 30 points grid in momentum space suggest that this does not spoil the 
reliability of the size of the effects we found. 

With reference to Table HI a few comments are in order. First, note that for sufficiently small values of £,e,^n all 
effect of sterile states on BBN is due to the increased A'cfi, which is in any case larger than what found in the single 
momentum approximation for the same value of This holds true for both Yp and ^H. On the other hand, for 
high values, say |Ce,^| ~ 10~^, some fraction of the sterile neutrino population builds-up relatively late, namely after 
the freeze-out of the active neutrinos. These cases are associated to a AiVeff significantly smaller than 1. While this 
quantity is still the main cause of the change in ^H, a significant fraction of the effect on Yp is due to the changes 
of the weak rates regulating the n p chemical equilibrium due to distorted i^e and distributions. In the case 
Ce = = 10~^, this effect is ~ 75% than the one related to the speed-up of the expansion due to AA'cfi, while 
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FIG. 3: Cases with = 10 ^ . Upper panels: Final Ve energy spectra at T = 1 MeV (dashed curve) and initial ones (continuous 
curve). Lower panel: Ratio R between final and initial i^e energy spectra. Left panels refer to = C/j while right panels are 
for = -Cm- 



for — = 10~^ it becomes three times larger than the other ^ . Phenomenologically, these scenarios with large 
asymmetries in the presence of sterile neutrinos yield comparatively lower values of A'eff as probed by CMB (say, 
iVoff — 3.2), while altering Yp by an amount loosely equivalent to a larger TVcff. We note that in these cases the 
presence of eV sterile neutrinos pushes the change of Yp in the direction of an increase with respect to the standard 
BBN value. This behavior confirms the prediction based on analytical estimations presented in [70[. The comparison 
with the last three rows in Table U also shows that this trend is opposite of what obtained in presence of a positive-sign 
electron neutrino asymmetry (in absence of sterile neutrinos). Actually, in some previous phenomenological analyses, 
not treating sterile neutrinos dynamically (like in (2^) it has been envisaged the possibility "to mask" the presence of 
extra eV-scale sterile neutrino degrees of freedom to BBN by introducing large chemical potentials, treating the system 
as a "degenerate BBN" plus extra radiation. This prescription may lead to even qualitatively wrong conclusions, such 



as the positive correlation in allowed regions between the increase of ^ and Ncr, visible e.g. in Fig. 6 of Ref. [22 1. 

Quantitatively, we observe that the modifications of BBN yields induced by sterile neutrinos are sizable, as reported 
in Table I. For comparison, the statistical error on the astrophysical Yp determination can be as small as 0.001 (albeit 
the systematic error is at the moment several times larger) [7l[, the determination of ^H/H in the highest quality 
Quasar system is (2.535 ± 0.05) x 10"^ f23|, and the 1 (T errors on N^s and Yp reported by the combined analysis of 
the Planck team amount to AA^off — 0.27 and AYp ~ 0.021 [26], in substantial agreements with earlier forecasts [t^ . 

Finally, we checked that changing squared mass differences and mixing angles within current uncertainties [65j , one 
can easily obtain 0(10%) differences in A^eff- Also, the results presented here are rather on the conservative side, as 
far as the particle mass parameter is chosen. For example, considering only disappearance experiments, larger values 
of sterile mass splitting than the 0.89 eV^ used here are preferred, see e.g. In this case, the results we obtained 
are modified in the sense of an easier thermalization of the sterile state. 



^ The evolution for this latter case is extremely slow and there is still a small evolution in the parameters taking place at T < 1 MeV. 
Hence the results presented here, which assume that the asymptotic results are equal to the ones at the smallest temperatures followed 
are slightly conservative. The actual effect should be a bit larger. 
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FIG. 4: Evolution of AA^cff vs. temperature T for equal asymmetries Ce = Cm C^ft panel) and opposite asymmetries Ce = —Cm 
(right panel). Solid curves refer to Ce = 10~^ while dashed curves are for Ce = 10~^. 



Case 


AiVeff 


cri 




^H/H (xlO^) 


Id < 10-' 


1.0 


1.0 


0.259 


2.90 


Ce = -Cm = 10-' 


0.98 


0.89 


0.257 


2.87 


Ce = Cm = 10-' 


0.77 


0.51 


0.256 


2.81 


Ce = -Cm = 10-' 


0.52 


0.44 


0.255 


2.74 


Ce = Cm = 10-' 


0.22 


0.04 


0.251 


2.64 


Ce = ICmI = 10-', no Us 


~ 




0.246 


2.56 


Ce = ICmI = 10-', no Us 


~ 




0.244 


2.55 


standard BBN 








0.247 


2.56 



TABLE I: The values of ATVcff and the calculated abundances of ^He mass fraction Yp and deuterium in the different cases 
considered in this paper. For comparison, the third column refers to the increase in the effective degrees of freedom obtained 
in the average momentum approximation, AN^^ . 

IV. CONCLUSIONS 

In this article, we have presented for the first time the results of a multi-flavour, multi-momentum computation 
of cosmological neutrino spectra in presence of a sterile state wi th p arameters consistent with those invoked in the 
interpretations of short-baseline neutrino oscillation experiments [1^ . We have considered the case of relatively large 
neutrino asymmetries {\Li,\ ^ 10"'^ — 10~^), which suppress the sterile production in the early universe giving in 
this way a better agreement with several cosmological observations. Our results show an enhancement in the sterile 
neutrino abundance of up to 0.2 effective degrees of freedom with respect to what observed in the average momentum 
study. This shifts the asymmetries needed for a significant suppression of N^s to relatively larger values, of the 
order of \L^\ ^ 10~^. Moreover, starting with opposite asymmetries in different fiavours with vanishing net neutrino 
asymmetries provides a less effective inhibition. On the other hand, for large asymmetries the significant production 
of the sterile neutrinos after the Ve and Pe decoupling causes non-negligible distortions (from few % up to ~ 20% 
for some cases considered) in the Ve and energy distributions. As a result — leaving apart the problem of how to 
generate such large initial asymmetries — the modifications of BBN yields are sizable, as reported in Table HI Also 
worth noting, while for ^H/H these modifications are essentially due to the larger NcS, for the "^He mass fraction Yp a 
significant effect follows from (anti) neutrino distribution distortions produced in presence of both asymmetries and 
a sterile state. 

Although it is not the purpose of the present paper to provide cosmological constraints on sterile neutrinos, we can 
outline the two following possible scenarios: 

i) thermalization of the sterile state in the presence of small or vanishing neutrino asymmetries. In the "interest- 
ing" parameter space this is often close to total, and already with pre-Planck data tensions were manifest between 
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cosmological disfavored regions and laboratory claims (see e.g. Fig. 5 in [t^ or the left panel of Fig. 3 in [s^l). This 
is amenable to a systematic analysis of the sterile neutrino mass and mixing parameter space, and a first analysis 
following the Planck data release is strengthening the above-mentioned tensions [t^ . 

ii) No complete thermalization of sterile neutrinos. Extra ingredients (which imply major cosmological changes) 
must be introduced. The most popular proposal in the literature has been to invoke neutrino asymmetries, whose 
effects on cosmological observables has been treated with simplified recipes. Here we put this possibility under 
close scrutiny, finding notable differences with naive expectations found in the literature. This scenario could lead 
to a possible inconsistency in the value of A^eff extracted from CMB and BBN. Indeed, sterile neutrinos and large 
asymmetries would produce a relatively low value of N^s as probed by CMB, and an increase of Yp. This latter would 
be mimicked by a larger A^off- It is also worth stressing that the standard BBN prediction Yp = Yp{uji,, N^tf) is not 
valid anymore. Given our current CPU power, we have limited our analysis of this scenario to few representative 
cases. Contrarily to what proposed in the scenario i) a large scan over parameter space (which would be needed for 
a quantitative statistical analysis combining CMB, LSS and BBN data) appears a hard task. Indeed, any precise 
quantitative result for the cosmological observables would depend on an interplay of active-active and active-sterile 
neutrino oscillations (multi-fiavour effects), whose matrix structure is not fully determined by current data on short- 
baseline neutrino anomalies (most notably in the r — s sector) . Hence "precision" computations would be illusory and 
premature, given the dependence from unknown or poorly constrained parameters. 

Future directions of these studies would depend on the fate of the short-baseline neutrino anomalies and on the 
current generation of cosmological measurements. Definitely, these results are already orienting — and will orient even 
more in the near future — further research on the non-trivial role played by light sterile neutrinos in cosmology. 
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